Scalable continuous- variable entanglement of light beams produced by optical 

parametric oscillators 
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We show that scalable multipartite entanglement among light fields may be generated by optical 
parametric oscillators (OPO). The tripartite entanglement existent among the three bright beams 
produced by a single OPO - pump, signal, and idler - is scalable to a system of many OPOs by 
pumping them in cascade with the same optical field. This latter serves as an entanglement distrib- 
utor. The special case of two OPOs is studied, as it is shown that the resulting five bright beams 
share genuine multipartite entanglement. In addition, the structure of entanglement distribution 
among the fields can be manipulated to some degree by tuning the incident pump power. The scal- 
ability to many fields is straightforward, allowing an alternative implementation of a multipartite 
quantum information network with continuous variables. 

PACS numbers: 03.67.Mn, 42.50.Dv, 42.65.Yj 



I. INTRODUCTION 

Entanglement is the primary resource in the field of 
quantum information. Its strikingly singular properties, 
which permeated the early debates on the fundamentals 
of quantum mechanics [T] , are nowadays employed to take 
computation technology and communication security be- 
yond the classical limit 0, Some applications have 
been successfully demonstrated in systems composed of 
up to four parties 0, [^, but their unrestricted use in 
quantum computation and networks relies on the ability 
to entangle numerous physical systems @. 

Multipartite entanglement brings a new level of com- 
plexity, and therefore novel possibilities, such as stronger 
violation of Bell-type inequalities The geometry of 
entanglement distribution among parties can be manipu- 
lated to perform one-way quantum computation [1, [9| or 
to guarantee that confidential information is correctly ad- 
dressed and safely transmitted [l^, [ll| . When it comes to 
continuous- variable (CV) entanglement [l^ , the quadra- 
tures of light fields - equivalent to the intensity and phase 
in the case of bright beams - offer the possibility to im- 
plement such ideas. 

The usual way to produce multipartite CV entangled 
light beams relies on the interference of at least one quan- 
tum field, presenting quadrature squeezing, with any de- 
sired number of vacuum fields using beam splitters |4l.ll3l|. 
This has been realized for up to four squeezed beams by 
employing optical parametric amplifiers [14| . However, 
the requirement for interference limits this type of scheme 
to fields with the same optical frequency. 

We present an alternative scalable way to directly gen- 
erate CV multipartite entangled beams without the need 
for interference [l^, [l^. Our scheme utilizes as a ba- 
sic resource the tripartite entanglement existent between 
the three fields - signal, idler, and pump - produced by 
an optical parametric oscillator (OPO) operating above 
threshold |17|. Various OPOs are used in a chain con- 



figuration, such that the pump beam refiected by the 
first one, already entangled to the first pair of down- 
converted twin beams, pumps a second OPO, entangling 
both pairs of twins, and so on. The pump beam then 
acts as a sequential entanglement distributor. Further- 
more, we show that the entanglement structure can be 
manipulated by tuning the pump power. We study the 
simplest case of two OPOs to show that the resulting five 
light beams present genuine pentapartite entanglement. 

This paper is organized as follows. We begin in Sec. lU 
with the method we use to characterize multipartite en- 
tanglement, the positiyity of the partially transposed 
(PPT) density matrix [S, We proceed in Sec. |TlT] 
and IIVI with a brief review of the equations describing 
the quantum properties of the light beams generated by 
a single OPO. By using this result, it is straightforward 
to derive the set of equations for the pentapartite sys- 
tem. In Sec. |Vl we present our results by applying the 
PPT criterion to the composite system of two OPOs, 
thus demonstrating the existence of genuine pentapartite 
entanglement among the light beams, and analyze the 
structure of entanglement distribution. Concluding re- 
marks are presented in Sec. IVII 



II. 



MULTIPARTITE ENTANGLEMENT OF 
GAUSSIAN STATES 



Gaussian states are completely characterized by their 
second order moments, organized in the covariance ma- 
trix V = (x^ x) , where 
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is the vector of the amplitude [pj = exp(— i(^j)aj -t- 
exp(z(/?j) a]] and phase [cjj — —i{ejip{—iipj)dj — 
exp(z(/?j) a])] quadrature operators, chosen relative to a 
phase ifj, and N is the number of field modes. The op- 
erators d, and d* are the usual annihilation and creation 
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operators for mode j. The canonical commutation rela- 
tions can be written in the compact form [x, jF] — 2if2, 
where 

n = ^u., co=(^^ M . (2) 
j=l ^ ^ 

A typical unitary physical transformation acting on 
the Hilbert space of the system, such as the ones as- 
sociated with beam splitters and one- or two-mode 
squeezers, preserves the gaussian character of the state. 
The corresponding transformation S acting on the co- 
variance matrix pertains to the real symplectic group 
S E Sp{2N,Ti), which preserves the commutation rela- 
tions, sns"^ = n^. 

In order to represent a physical state, the covariance 
matrix must respect the Robertson-Schrodinger uncer- 
tainty principle [2l|, 

V + in>0. (3) 

This condition, which is invariant under symplectic trans- 
formations, implies a constraint to the symplectic eigen- 
values of the covariance matrix, 

i^,>l, kel,2,...,N. (4) 

The symplectic eigenvalues can be computed as the 
square roots of the ordinary eigenvalues of —{flVy. 
They are invariant under local operations and classical 
communication (LOCC) [l^l- If the covariance matrix 
corresponds to a pure state, then one necessarily has 
i/fc = 1, Vfc. 

The existence of entanglement among the field modes 
may be unveiled by applying separability tests to all pos- 
sible partitions, and its structure, also by the reduced 
subsystems. If all partitions of a pure state are insep- 
arable, the system is said to present genuine TV-partite 
entanglement. One powerful testing method relies on the 
positivity of the partially transposed (PPT) density ma- 
trix [3]. This map is positive regarding all separable 
states, but may be negative to entangled states. It is 
equivalent, in phase space, to a local inversion of time 
for the transposed subsystems [l^. Regarding the co- 
variance matrix, the PT operation takes in for 
the desired subset of modes. The failure of the resulting 
PT covariance matrix V to comply with the uncertainty 
principle ([3]) is a sufficient condition for the existence of 
entanglement between the transposed subset and the re- 
maining subset [l^ . 

Therefore, Eq. [3] must be satisfied by all the possible 
transpositions of covariance matrices representing com- 
pletely separable states. Equivalently, the symplectic 
eigenvalues i> oi V must fulfill Eq. [?] in this case. The 
PPT criterion is both necessary and sufficient for pur e 
or mixed states split in partitions 1 x (A^ — 1) [l9l. l23l|. 
Other partitions from systems with N > 2 may possess 
bound entanglement, a non-distillable form of entangle- 
ment which is not revealed by partial transposition [23| |. 
Nevertheless, it is always sufficient. 



It is worth mentioning that the smallest symplectic 
eigenvalue i>„^i„ of the PT covariance matrix is useful not 
only to witness the entanglement but also to quantify it. 
In fact, the entanglement measure given by the logarith- 
mic negativity 24] can be written as a decreasing func- 
tion of i>,„i„, for all A^)-mode bisymmetric Gaussian 
states [22,]. Thus, a larger violation of Eq. (jH) implies a 
larger amount of entanglement. 

III. DESCRIPTION OF THE OPO 

The optical parametric oscillator consists of a non- 
linear x'^' crystal disposed inside an optical cavity, in 
this manner coupling three modes of the electromagnetic 
field. It is driven by an incident laser, the pump beam, 
at frequency lOq. Following the usual terminology, the 
fields generated by downconversion are called "signal" 
and "idler," with frequencies loi and respecting, by 
energy conservation, uJq = + uj2- Above the oscilla- 
tion threshold, signal and idler are bright light beams 
known to possess strorig quantum correlations in both in- 
tensity [25| and phase [2^ , therefore called "twin beams" . 

The covariance matrix of a single OPO can be found 
in many references (27l |28| . The standard treatment be- 
gins with the master equation for the density operator, 
which is then converted to a Fokker-Planck equation for 
a quasi-probability distribution (in the present case, the 
Wigner function), and finally to a set of Langevin equa- 
tions for its complex arguments representing the classical 
field amplitudes aj{t). The labels j S {0,1,2} corre- 
spond to pump, signal and idler, respectively. The set of 
Langevin equations describing the quantum fluctuations 
da^ of the intracavity fields {t) are: 

T^(5ao = -7o(l - «Ao)(5ao + y^oSa"^ - 

-—7(1 + zA)e-'''(e'''2(5ai + e'-^^Sa^), (5) 

Po 

T—Sa^ = —7(1 — iA)Sai + \/2r/Svi + 
at 

■yd — iA) 

+ — -e'''{pe-"'^Sao+Poe'''°Sa;), (6) 

Po 

T—5a2 -7(1 - iA)5a2 + ^J^qbv^ + 

-f^^^^— ^e"'(pe-"'Mao-hpoe'"'''fc;), (7) 

where a linearization procedure has been applied, ckj (t) = 
aj-^6aj(t\ such that terms involving the product of fluc- 
tuations have been ignored. The mean complex ampli- 
tude (5tj = exp(i(/9j) reproduces the classical result for 
the mean intracavity field amplitude p^ and phase psi ] 
(in the above equations, Q = i-p^ + '-Pi — fo)- The latter 
is defined with respect to the incoming pump beam aj,", 
chosen real. The coefficient 7^ is related to the coupling 
mirror transmission for the field j, — 2^^ <C 1. Only 
intracavity losses originating from mirror transmissions 
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are considered in our analysis (although the effects of 
spurious losses are briefly discussed in the final results). 
The detuning between the OPO cavity and the field 
is given by Aj, where Ai = Aa = A, and the cavity 
roundtrip time is t. We assume equal mirror transmis- 
sions for signal and idler beams, i.e. 71 = 72 = 7, result- 
ing in the equality of their mean amplitudes, = p2 = p. 
The terms 6vj are the vacuum fluctuations associated to 
these couplings. For the intracavity pump mode, they are 
due to the quantum fluctuations of the input pump laser 
beam, Svg = JaJ,"; however, since we consider it as shot- 
noise limited, its fluctuations obey the same statistics as 
the vacuum. 

To obtain from Eqs. ©-([T]) equations for the quadra- 
tures, we write 



SpM+^SqM] 



(8) 



where the phase references for the definition of the 
quadratures were chosen as the respective mean fields' 
phases. We further simplify the problem by assuming 
exact triple resonance, Aq = A = 0, since detunings 
have the main effect of weakly coupling amplitude and 
phase quadratures, but do not change the overall physical 
behavior of the system 28]. The solution is readily ob- 
tained in frequency domain by considering the combined 
quadratures p± = {pi ± p^l\/2 and = (gi ± q.})l\pl, 
for which the evolution equations decouple [27| . The out- 
put quadratures are finally determined using the input- 
output relations. 



(9) 



6pT{^) = -K.(f^) + ^J'2.-ibpX^), 



where s € {0,+,—}. A similar relation holds for the 
phase quadrature. The solution is 



SpT 

SpT 

Sq-y 



(5pi," -I- (27^p - 1) 5v^^ 
Sq^" + (27^, - 1) (5w,^ 

I?, Sql" - 6vg^ , 
in' ^ 

-6v„ , 



l + ifl' 



where 
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27^/3^ 



7o + 2i-fil' 



^p 



7o - 
-1 + 



2^717' 
270 



7o -I- 2i^n' 



1 



(10) 

(11) 
(12) 
(13) 

(14) 
(15) 

(16) 
(17) 

(18) 
(19) 



The parameter fJ' = n/duj is the analysis frequency 
relative to the OPOs cavity bandwidth for the twin 



beams, and /? = p/po- The expressions for 5pi 2 and 
Sqi 2 are directly obtained from the appropriate lin- 
ear combinations of the above equations. The covari- 
ance matrix is calculated as ^ = {x^ x), with x ~ 
{5pi , 5qi , 5p2 , 5q2 , 5po , Sq^) (the superscript "out" has 
been dropped for the sake of notational simplicity). 




Vacuum 



FIG. 1: (Color online) Elementary symplectic transforma- 
tions bringing a three-mode vacuum state into the fields pro- 
duced by an OPO: pump (0), signal (1) and idler (2). A beam 
splitter transformation is represented by BS and a single- 
mode squeezer by Z. The circles and ellipses represent the 
fields' quadrature noise powers in phase space. The existence 
of tripartite entanglement among the beams 0, 1, and 2 is 
clearly seen as arising from the interference of squeezed beams 
in this equivalent scheme. 

The main characteristics of this three-beam system 
were shown to be the following. The twin beams are 
entangled due to strong intensity correlation and phase 
anti-correlation (25l . [2g| . Because of pump depletion 
which always occurs above threshold, the twin beams are 
able to influence back the pump field, occasioning entan- 
glement between the sum of twins and the reflected pump 
field. As a consequence, the OPO directly produces tri- 
partite entangled light fields [13, [2§| . 



IV. TRIPARTITE ENTANGLEMENT 

The three beams generated by a single OPO form a 
pure quantum system, and hence can be described by 
a symplectic transformation acting on the three-mode 
vacuum field, Vs = S'opol<S'Jp„, where 1 is the 6x6 identity 
matrix. In order to provide an intuitive picture of the 
tripartite entanglement. Fig. [T] illustrates that the above- 
threshold optical parametric oscillator can be thought of 
as a device which realizes the symplectic transformation 
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(20) 



(21) 



Z,,,, =diag[e"i',e ''>',e''j',e "3'], (22) 

are the matrices describing the beam sphtter and the 
squeezer transformations, respectively. The subscripts 
i',j' S {0, 1, 2, +, — } indicate the subspace to which the 
transformation is apphed; absence of subscript means 
that it acts on the complete system. The squeezing pa- 
rameter r is such that the same amplitude and phase 
noise power occur in both modes. In addition, the OPO 
displaces the fields in phase space by large amounts a^; 
however, these operations affect only the classical mean 
amplitudes, and are therefore ignored in the treatment 
of the fields' quantum properties. 

The existence of entanglement and its structure (i.e. 
the way it is distributed among the fields) is revealed by 
the PPT criterion. We employ Eqs. pi?)) -- p5|) to con- 
struct Vs • The partial transposition operation may then 
be applied to either the pump beam, resulting in the 
PT covariance matrix Vi"', or to one of the twin beams, 
yielding Vg'^' = ^i^' tilde indicates partial transposi- 
tion with respect to the field mode labelled by the super- 
script). For this type of partition, at maximum one sym- 
plectic eigenvalue vlf^ from each PT matrix may assume 
a value smaller than one [30| , in this case demonstrating 
entanglement by violation of Eq. Q. Therefore, we al- 
ways consider the smallest symplectic eigenvalue i''^' to 
witness entanglement. 

We employ in our numerical results the typical experi- 
mental values 7o = 0.05 and 7 = 0.01 for the OPO cavity 
mirrors transmissions, analysis frequency Q' — 0.5 (rela- 
tive to cavity bandwidth) and pump power cr = 1.5 (rela- 
tive to threshold) . The entanglement between pump and 
twins is quantified by = 0.43 < 1. In the same man- 
ner, each twin beam is found to be strongly entangled to 
the remaining two beams, since = — 0.20. These 
two eigenvalues are very small close to threshold (u « 1) 
and remain smaller than one even for high pump power. 
On the other hand, i>'-°^ reaches one at threshold, having 
its minimum (maximum entanglement) around the cho- 
sen value of cr = 1.5 [l3|- This is physically explained by 
the fact that very close to threshold all quantum correla- 
tions reside between the twin beams; on the other hand, 
at cr = 1.5 the three output fields exhibit similar output 
mean power, usually the best situation for the equal shar- 
ing of correlations. Finally, all three symplectic eigenval- 
ues are very small (i' — ^ 0) for null analysis frequency 
and approach unity as Q' increases, since pump-signal- 
idler quantum correlations are well known to reach their 
maxima at lower frequencies. These results are enough to 
characterize the tripartite entanglement as genuine [sij . 
but the structure of entanglement is further revealed by 
the reduced covariance matrices. 

Tracing out the pump beam, we find that the PT re- 
duced matrix V12 (signal and idler) still possesses a small 
symplectic eigenvalue, i>i2 — 0.28 < 1. Thus, signal 
(idler) is strongly entangled to the remaining two beams, 
as shown before, because it is actually highly entangled to 
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FIG. 2: (Color online) Tripartite entanglement structure 
for pump, signal, and idler (circle, up-triangle, and down- 
triangle, respectively). The numbers represent the smallest 
symplectic eigenvalues of the PT matrix connecting the in- 
dicated subsystems. On the left side, entanglement between 
each pair of beams tracing out the third is presented. The 
overall structure of entanglement among the three beams is 
illustrated on the right, with pump entangled to the signal 
and idler highly entangled system. 

idler (signal). On the other hand, the PT reduced matrix 
for the bipartite pump and idler or pump and signal sys- 
tems, V(n = V02, possesses a symplectic eigenvalue close 
to 1, Doi = ^0-2 = 0.89, suggesting a reduced degree of en- 
tanglement between pump and a single twin. In fact, for 
a broad range of experimental parameters, the system 
presents Uoj ^ Ij but, contrarily to the other symplec- 
tic eigenvalues, it tends to one as the analysis frequency 
r2' is set close to zero. This result is again understood 
in terms of the increasing correlations between twins as 
r2 — > 0, such that 0, implying perfect entanglement 

between twins beams and no entanglement with pump in 
this limit. 

It follows from this discussion that the pump beam is 
more entangled to the signal and idler composite subsys- 
tem than to either signal or idler alone (i.e. < zJqi)- 
One concludes that the tripartite entanglement in the 
OPO is not just a consequence of three pairs of bipartite 
entangled subsystems, since there is a finite amount of en- 
tanglement recoverable only in the complete three-beam 
system. Fig. [5] pictorially represents the distribution of 
entanglement among the pump, signal, and idler beams 
for a single OPO. 



V. SCALABLE CV ENTANGLEMENT 

We now consider a second OPO (designated "OPO i?" ) 
pumped by the beam reflected from the single OPO of 
last section ("OPO A"), in a chain configuration (Fig. [3]). 
The pump power reflected by the flrst OPO (-Pq""*'*) is a 
decreasing function of the incident pump power (P"), 
since a larger fraction of its energy is down-converted 
for increasing input power. Therefore, the second OPO 
must have a lower threshold power (P^f ) than the first 
(Pj'jJ) in order to operate above its oscillation thresh- 
old. This could be accomplished either by choosing ap- 
propriate mirror transmissions or nonlincarity strength 
X^, since Pti, « 7o7^/x^- It turns out that the ra- 
tio between the thresholds is the most relevant factor 
for the entanglement between the two pairs of twins, 
as we will show later. Therefore, as a a = P'"/Ptt in- 
creases, (Tb = Po°"''*/Pth decreases, following the relation 
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Reflected pump 
(label 0) 



FIG. 3: (Color online) Sketch of the system considered. Two 
OPOs are pumped in a cascaded configuration, generating 
scalable entanglement among five beams. Pump beam dis- 
tributes the entanglement. 

The pentapartite system is composed of the two pairs 
of twin beams generated by the two OPOs (denoted by 
1^, 2^, 1b, and 2b), and the pump beam reflected by the 
second OPO (denoted by 0). We define the covariance 
matrix Vc — {x^ Xc) for the composite system, with 



{5pt Sqf, 5pi 
Spf, Sqf, 6p^, 



hi 



(23) 



Spo, Sq„). 



The twin fields produced by the OPO A respect ex- 
actly the same solutions of Eqs. (fT0 l) -(fT5 l) . such that 
{5pf,Sqf} = {dp,,Sq^,}. The solution for the OPO B 
is found in a straightforward way. As the subtraction of 
twin quadratures does not depend on the input pump, 
6p^ and Sq'^ are also given by Eqs. (fT4|) .(|15 p . As for the 
sum and refiected pump subspaces, the same solutions 
of Eqs. lfT0)) - p3l) apply for Spf and Sqf using as input 
pump instead the reflected pump from the OPO A, i.e. 
by the substitutions {(5po",(5g"} {Sp^,Sq^}. In this 
manner, quantum correlations from the pump beam re- 
flected from the OPO A are transferred to the twin beams 
produced by the OPO B. 

In order to demonstrate genuine pentapartite entangle- 
ment, it is sufficient to show that all possible bipartitions 
of the system are entangled, since the complete system 
is pure. We use the numerical values 7^ = 0.04 and 
7^ = 0.0075 for the cavity mirrors of OPO B (the param- 
eters for OPO A are the same used in the last section), 
and P^/P^ — 0.45. We begin by studying the system 
behavior as a function of the incident pump power at the 
fixed analysis frequency U,' — Q.\ (relative to the cavity 
bandwidth of OPO A). 

First of all, we investigate how each field is entangled 
to the remaining ones (partitions of the form 1x4). The 
set of the smallest symplectic eigenvalues v^^^ resulting 
from the partial transposition with respect to each 
one of the five beams presents the following properties. 
The first result is that each of the four down-converted 
beams is highly entangled to the remaining four beams, 
j)'.^A)^i)(2A)^i)(iB)^i}(^B) ^ I This was to be expected, 
since strong entanglement is actually present inside each 
pair of twin beams. The interesting result regards the 
pump beam, for which t'*"' <C 1 (Fig. 21 full squares) 



attests that its strong entanglement to the two pairs of 
twins remains. 
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FIG. 4: (Color online) Smallest symplectic eigenvalues ob- 
tained by applying the PT operation to various partitions of 
the total covariance matrix. The superscripts indicate the 
transposed subsystems, while subscripts, when present, indi- 
cate which reduced system is being considered. A value below 
one proves the entanglement between the subset of transposed 
subsystems and the remaining ones. 

In the same manner, we consider partitions of the type 
2x3. Fig. [4] summarizes the results. Entanglement be- 
tween each pair of twins and the remaining three beams is 
attested by the smallest symplectic eigenvalues obtained 
by transposing either subspace A or B, < 1 (full cir- 
cles) or < 1 (open circles). One sees that i>'-"^ 
when ffg ~* I, and similarly u^"'' —^ i)'^' when cr^ — )■ 1, 
as expected for the limit case when just one OPO oscil- 
lates. We note that the value of £'*'^' at O' = is the 
same as previously obtained by transposing the reflected 
pump beam from a single OPO (denoted as in last 
section): the twins from OPO A are as entangled to the 
remaining subsystem of three beams as they were to the 
single OPO reflected pump in Sec. IIVI indicating that 
the original entanglement between twins A and reflected 
pump A is redistributed to the three new beams. This 
intuitive picture is not strictly valid anymore for higher 
analysis frequencies, for which the original entanglement 
can decrease slightly. Finally, we have observed that the 
symplectic eigenvalues slightly decrease for lower values 
of Q' . In addition, their behaviors do not depend signifi- 
cantly on the other experimental parameters, implying a 
robust pentapartite entanglement. 

These results suffice to prove the existence of penta- 
partite entanglement in the system. It is not necessary 
to test partitions of the kind "signal A and pump related 
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FIG. 5: (Color online) Entanglement distribution among the 
five beams: twins A (open triangles), twins B (full trian- 
gles), and reflected pump (circle). The full arrows connect 
subgroups which violate the PPT criterion, while the dashed 
ones point out that the entanglement structure changes with 
n', a, and PfjP.t 



to the remaining three beams," since we know signal A 
is highly entangled to idler A and, as a consequence, to 
the system of remaining beams. The symmetry included 
in the system by the existence of two pairs of highly en- 
tangled beams reduces the number of effective partitions 
which need to be tested. The important results for now 
are the entanglement between pump and the remaining 
four beams, and between each pair of twins and the three 
remaining beams. 

The entanglement structure is again better revealed by 
the reduced covariance matrices. The results are sketched 
in Fig. m We trace out pump subspace to obtain the re- 
duced covariance matrix Vab for both pairs of twins, and 
apply the PT operation on one pair of twins. In Fig. 21 
the smallest symplectic eigenvalue i>'j^ < 1 (open trian- 
gles) attests the entanglement between the two pairs of 
twins for cr^-i < 1.65. Therefore, the pump beam reflected 
by the first OPO can effectively entangle both OPOs out- 
puts. The eigenvalue tends to one from below (i^Ib' 1) 
as a A — > 1; for increasing pump power, it changes from 
i^Aji < 1 to j^ie' = 1) modifying the entanglement struc- 
ture. 

The next reduced subsystem Vgo excludes the twins A 
to probe the entanglement present between pump and 
the twins B. The partial transposition of the pump field, 
resulting in z/bo' < 1 (open squares), shows that they 
are always entangled. Its value increases (weaker en- 
tanglement) as the threshold of OPO B is approached 
(ctb 1), which is to be expected for a single OPO 
(SecHVl). 

Finally, the reduced covariance matrix Vao of pump 
and twins A is transposed relative to pump in order to 
probe their entanglement. It turns out that Vao fulfills 
the uncertainty principle for lower input pump powers 
{a A < 1.6), since i^ao = 1 (full triangles), meaning that 
pump is not entangled to the twins A alone in this region. 
In view of the previous results, this indicates that the en- 
tanglement between twins A and the pump beam prior 
to interaction with OPO B is strongly converted into en- 
tanglement between the pairs of twins. At the same time, 
the scalability is warranted by the entanglement appear- 
ing between twins B and the final reflected pump. On the 
other hand, for cr^ > 1.6 the pump remains entangled to 




■f fli^ ■ 'Ptii'^ (ratio between tlireshold powers) 

FIG. 6: (Color online) Contour plot for the smallest symplec- 
tic eigenvalue considering the reduced subsystems: pump and 
twins A (patterned region) and twins A and twins B (region 
without pattern). The structure of pentapartite entanglement 
changes according to Q' and P^. The small black region cor- 
respond to a situation where both u')^ and u^^q are below 
one, although with a weak entanglement. It has been used 
input pump power a a = 1.1. 



twins A (as for a single OPO), but fails to transfer these 
correlations to the two pairs of twins, which remain only 
indirectly connected to each other by the entanglement 
they share with the pump beam. 

Thus, as the incident pump power is varied, an alter- 
ation in the entanglement structure takes place, switch- 
ing between D"}^ < 1 and D^o < 1- One may choose 
appropriate parameters (analysis frequency il' and OPO 
B's threshold power) in order to have one or both of the 
entanglement types: between pump and twins A or be- 
tween the two pairs of twins. Fig. [S] depicts this behavior 
at the flxed input pump power a a = 1.1. The region 
without pattern corresponds to j^Ib' < 1 and the pat- 
terned region to ufo < 1- Stronger entanglement (darker 
tones) between the pair of twins is found for small val- 
ues of Q' and P^'^/Pit: while pump and twins A get more 
entangled as the OPOs become similar (P^f « Ptt)- 
small coexistence region, although presenting weak en- 
tanglement (1 > C'il' « jyfd ^ 0.90, seen at tr^ « 1.63 
in Fig. 2]), is depicted by the black area. We also note 
that it is possible to adjust the ratio between the OPOs 
threshold powers either by selecting adequate crystals or 
cavity mirrors, since these results are not much sensitive 
to the exact choice of mirror transmissions. 

As a last remark, we mention that the inclusion of 
small spurious losses (< 0.05 x 70,7) does not change 
the qualitative behavior of the system. The presented 
symplectic eigenvalues become larger, indicating a loss 
of quantum correlations among the beams. Intracavity 
losses in the twins modes have an important quantitative 
effect in this direction, specially concerning the entan- 
glement between the pairs of twins. Spurious losses in 
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the pump mode have the main consequence of decreas- 
ing the available power for pumping the OPO B, effec- 
tively reducing the horizontal axis range of Fig. |4] (with- 
out altering its qualitative features). The most affected 
symplectic eigenvalue in this case is again related to the 
entanglement between the pairs of twins. In both cases, 
a limited region in a a then appears in which neither the 
reduced subsystems of pairs of twins nor twins A and 
pump are entangled. Nevertheless, full inseparability is 
still attested. 



VI. CONCLUSION 

The entanglement existent among the three fields pro- 
duced by an above-threshold optical parametric oscillator 
allows the implementation of a scalable network of multi- 
partite entangled light beams, by successively employing 
the pump beam reflected by one OPO to pump a sub- 
sequent one. The pump beam acts as an entanglement 
distributor among otherwise independent OPOs. 

Moreover, the internal structure of the final penta- 
partite entangled state can be manipulated to some ex- 
tent by tuning the incident pump power. The analysis 
of the reduced covariance matrices considering two cas- 
caded OPOs reveals situations where the pump beam is 



either entangled to each pair of twins considered alone or 
both pairs are entangled to each other. In any of these 
cases the system is scalable, since each beam is always 
entangled to the remaining four beams. 

The major limitation concerning the maximum num- 
ber of beams which can be entangled comes from the drop 
in pump power as more OPOs are added to the multi- 
partite system. However, by choosing high-quality op- 
tics (low spurious losses) and periodically poled crystals 
which allow a very low oscillation threshold, the number 
of entangled beams could be increased. 

One special feature, the involvement of nonlinearities, 
distinguishes this system from other CV multipartite 
state generation based on interference, since this scheme 
entangles various spectral regions and, therefore, allows 
quantum information to be conveyed among different 
parts of the spectrum. Multicolor quantum networks of- 
fer the possibility to communicate quantum hardwares 
with otherwise incompatible working frequencies. 
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